liz22_del_14_33

seascape <- seascape_probability(reefoutline=reef_map, habitat=benthic_map)


liz22_del_14_33 <- st_read("/Users/rof011/coralseed/code/2021_2022_HighPRT_sim_examples/run_day_12038_liz22_del_14_33/run_day_12038_liz22_del_14_33.json")
## Reading layer `run_day_12038_liz22_del_14_33' from data source 
##   `/Users/rof011/coralseed/code/2021_2022_HighPRT_sim_examples/run_day_12038_liz22_del_14_33/run_day_12038_liz22_del_14_33.json' 
##   using driver `GeoJSON'
## Simple feature collection with 605000 features and 3 fields
## Geometry type: POINT
## Dimension:     XYZ
## Bounding box:  xmin: 145.4446 ymin: -14.71132 xmax: 145.4806 ymax: -14.69039
## z_range:       zmin: 3.55 zmax: 3.55
## Geodetic CRS:  WGS 84
particles <- seed_particles(input=liz22_del_14_33, seascape=seascape, probability="additive",
                            limit_time=24, competency.function = "exponential", 
                            simulate.mortality = "typeI", simulate.mortality.n = 0.1, 
                            return.plot=TRUE, silent=FALSE)

## ## coralseed model summary: #### 
## Importing 5000 particle tracks
## [No seed set, random draws used] 
## 
## Time start = 2022-12-18
## Time end = 2022-12-19 00:57:00
## Total dispersaltime (hrs) = 1.03958333333333
##   
## Competency at t6 = 639 / 5000 larvae 
## Competency at t12 = 1182 / 5000 larvae 
## Competency at t24 = 2100 / 5000 larvae 
##   
## Survivorship curve  
## Mortality at t6 = 19 / 5000 larvae 
## Mortality at t12 = 120 / 5000 larvae
settlers <-  settle_particles(particles, probability="additive", silent=FALSE, return.plot = TRUE)
##   
## 2061 / 5000 larvae settled 
## 
map_coralseed(seed_particles=particles, settle_particles=settlers, 
              seascape_probability=seascape, subsample=500, restoration.plot=c(100,100), 
              show.tracks = TRUE)